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We study in details the space-time dependence of the production of muonic, pionic, and other 
exotic atoms by the coherent photon exchange between nuclei at the Large Hadron Collider (LHC) 
at CERN. We show that a multipole expansion of the electromagnetic interaction yields an useful 
insight of the bound-free production mechanism which has not been explored in the literature. 
Predictions for the spatial, temporal, and angular distribution, as well as the total cross sections, 
for the production of exotic atoms are also included. 
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I. INTRODUCTION 

It is undeniable that ultra-peripheral collisions (i.e., collisions dominated by the electromagnetic interaction) be- 
tween relativistic heavy ions allows interesting physics to be probed. As far back as 1989 it was noticed that even the 
Higgs boson could be produced at comparable rates as in central collisions between relativistic heavy ions [T, '2] , the 
advantage being that ultra-peripheral collisions are cleaner in the sense that nothing else than the Higgs boson would 
be observed. By now, ultra-peripheral collisions are known as and excellent tool for several interesting phenomena 
and have been discussed in numerous publications, reviews, and popular articles (see, e.g. 00]). 

A process of interest in ultra-peripheral collisions with relativistic heavy ions is the production of pairs in which 
one of the particles is captured in an orbit around one of the ions in the collider ( "bound- free" pair production). In 
particular the process of pair-production with capture was used as a tool to generate the first artificial anti-hydrogen 
atom in the laboratory ^J. The pioneer CERN experiment was followed by another experiment at Fermilab 3]- The 
number of anti-atoms produced was shown to be in good agreement with detailed theoretical predictions |7]. 

In this article we extend these studies by considering the production of exotic atoms in pp and heavy ion collisions 
such as those being carried out at the Large Hadron Collider (LHC) at CERN. In our notation, an exotic atom is an 
atom in which one electron has been replaced by other particles of the same charge. In contrast to previous approaches, 
we show that a multipole expansion of the electromagnetic field yields several insights in the production mechanism 
which have not been explored in the literature. Our calculations apply to the production of muonic atoms, pionic 
atoms, protonium, etc. The Bohr radius for a muonic atom is much closer to the nucleus than in an ordinary atom, and 
corrections due to quantum electrodynamics are important. The energy levels and transition rates from excited states 
to the ground state of muonic atoms also provide experimental tests of quantum electrodynamics. Hadronic atoms, 
such as pionic hydrogen and kaonic hydrogen, also provide interesting experimental probes of the theory of quantum 
chromodynamics. A protonium is antiprotonic hydrogen, a composite of a negatively charged antiproton paired with 
a positively charged proton or nucleus. Protonium has been studied theoretically mainly by using non-relativistic 
quantum mechanics, which yields predictions for its binding energy and lifetime. The lifetimes are predicted in the 
range of 0.1-10 microseconds. While protonium production has a very small cross section in pp collisions at CERN, 
the cross section is appreciable for heavy ions, e.g. for Pb-Pb collisions. In this case, the anti-proton will be captured 
in an orbit around one of the Pb nuclei. Such a system would be of large interest for understanding inter-nucleon 
forces in charge-conjugate channels. 

Most of our calculations will be for the production of muonic atoms and pp collisions. There is no qualitative 
difference for the production of other exotic atoms via the same mechanism, except for the obvious reduction of the 
production yields due to mass differences. The production cross sections for Pb-Pb collisions will be enhanced by a 
huge amount, approximately equal to lO^'^ — 10""^^, but the details of the production mechanism are practically the 
same as in pp collisions. In Section II we show how a multipole expansion allows for a clear space-time description 
of the bound-free production mechanism. Useful approximate formulas are derived. In Section III we present our 
numerical results and discuss the physics properties of our results. Our conclusions are presented in Section iV. 
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II. PRODUCTION OF EXOTIC ATOMS IN ION-ION COLLISIONS 

In the frame of referenee of one nucleus, the time-dependent electromagnetic field generated by the other nucleus 
is given by the Lienard-Wiechert potential 

TZe^ 

A^,{Y,t) = v^cj), with = (l,v), and (/)(r,i) = ^y^, (1) 

where 

v = {x,y,z), R = (a;, y,rz), li! = {bx,by,Vvt), 

and V is the relative velocity between the nuclei which we will take to lie along the z-axis. T is the relativistic Lorentz 
factor r = 1/(1 — u^)^/^ (otherwise explicitly stated, here we use the units h = c=\). 

In first-order perturbation theory, the pair-production amplitude at time t for a collision with impact parameter 



b — Y + is given in terms of the transition density p(r) and the current density j(r), as 

ai,t(p,&,t) = - r diV"*' / dVj^(r)A''(r,t'), (2) 



where a; = e -f m — /, with e equal to the energy of the free positive particle and / being the ionization energy of the 
negative captured particle, m is the rest mass of either particle. The transition current is given in terms of the Dirac 
matrices 7^, and the particle and anti-particle wave functions, i.e., jp(r) = e\E'_7^\E'-|-, where is the wavefunction 
for the captured negative particle and 5*+ that of the free positive particle. 

The multipole expansion of j^A^^ has been extensively discussed in details in Refs. [8l-[T0]. Replacing the Schrodinger 
by the Dirac currents in their results, one finds that j^(r)A^(r, i) = X^Trin ^7r;K(r, 6, <), where ir = E,M and I — 
0, 1, 2, • • • , and k = —I, • • • , / denote the multipolarities. For electric El (electric dipole) and E2 (electric quadrupole) 
multipolarities, 

V,.A.,.,. i) ^ (r) (1 - ,r> (.) { (3) 



T7 2 / q /in I (if « = ±2) 

yE2.(r, b, t) ^ ^. (r) (1 - 7o7.)r^y2. (r) ^+ (r) ' ^ ^' { T{T^ + l)bvt, (if k = ±1) (4) 

(^,2 + \ ^273 {2TW ~ b^) (if K = 0). 



In addition, for magnetic dipole (Ml) [TU], 

TZe^y^2^ ( ±b, (if K = ±l) 



(5) 



For the positive particle wavefunction we use a plane wave and a correction term to account for the dis- 
tortion due to the nuclear charge [7|. The wavefunction is given by 4*+ = W [v(p) exp(ip • r) -I- 4''^] , where 
J\f{e) — exp(7ra_|_/2)r(l -|- ia-|_), with a+ = Ze^m/p, with p being the proton momentum, v(p) the particle spinor and 
|A/'(e)| = 27ra+/[exp(27ra+) -I- 1]. The correction term is given by eq. (B5) of ref. (7j, which for our purposes can 
be written as 



Kir) = f^v(p) / cf,e--^ Y'^^:^'t 

27r^ J (p-q)^(r-p^) 



(6) 



with e = \/ + Tii^ ■ Parts of this integral can be done analytically. 

For the negative particle we use the distorted hydrogenic wavefunction [71 for capture to the ground state. 



1 fZ^"^ 



4'_(r) = 



1 _ 1/^ . V 



u(eo) exp (-Zr/a) , (7) 



where u is the negative particle spinor, eq is the energy of the exotic atom bound-state, and a = 1/me^ being the 
hydrogen Bohr radius. As noted in [7), the reason why we need to keep the corrections in the wave functions to first 
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order in Ze^ is because for e ^ m these corrections yield a term of the same order of magnitude as the plane-wave 
to the total cross section. The reason for this are the small distances which enter in the calculation of the integral 
in Eq. ([2]). The wavefunction corrections are essential to account for the short distance effects properly. Such effects 
are not as important in the case of bound-free production of heavy particle-antiparticle pairs, but we keep them for 
a more accurate description of the process. 

To gain insight into the pair production probability with capture and its distribution in space-time, we will compare 
our numerical results with simplified calculations which neglect the above mentioned short-distance corrections of 
the wavefunction of the captured negative particle as well as of the free positive particle. This is a reasonable 
approximation for small values of the free particle energy, i.e. for e ~ m, as we show later. This approximation allows 
us to obtain the coordinate integral in a\st analytically, by using 



Jtx \ a 



3/2 



5/2 



{i + xiY 

6{a/Z)xl 



" yi,(p) for A = l, 



(8) 



Y2AP) for A = 2. 



where Xp = ap/Z and p is the unit vector for the positive particle momentum direction with respect to a z-axis 
pointing along the beam direction. To calculate total cross sections we keep the wavefunction corrections, Eqs. ([6]) 
and ([7]). In this case, some of the coordinate integrals need to be perfomed numerically. 

The exotic atom production probability density, at time t and for a collision with impact parameter h, is obtained 
by the square of expression ([2]), i.e.. 



(9) 



spins 



where the sum over spins is performed with standard trace techniques. 

The production probability, at time t and for a collision with impact parameter 6, is obtained by integration over 
momentum, 



P{b,t) 



(2^ 



r{p,b,t). 



(10) 



If one neglects the derivative corrections in the particle and anti-particle wavefunction, the sum over spins in eq. ^ 
yields 



_fp_ 
(27r)3 



y?T^A/'(e)FA«(p) 



|GAK(6,^,t)r, 



(11) 



where the small binding energy / was neglected as it is small compared to the mass m. The functions Gak are defined 
as Gak — Jq exp{—iujt'}{- ■ ■ }, where the terms inside braces {• • • } are the time-dependent terms to the right of the 
positive particle wavefunctions in eqs. (|3][5|. For the Ml magnetic multipole the same equation can be used with the 
replacement + m — >■ ^e — and multiply it by an extra factor {pz/mY . 

The cross section for bound-free pair production is obtained by integrating the production probability over all 
impact parameters at t = 00, i.e.. 



(fibP[b, 00). 



(12) 



The above integral over impact parameters diverge if the potentials of eqs. (|3])-([5]) are used. The reason is that these 
potentials are obtained from an expansion of the full Lienard-Wiechart potential of eq. ([T]) which is only valid for 
distant collisions, i.e., when the Lorentz-modified projectile coordinate R' is larger than the Lorentz- modified internal 
coordinate R. A better approach is to use a full representation of the potential in terms of a momentum transform, 
as done in ref. [7], 



I.e., 



, , TZe 
27r"^ 



■iq- [R-R'(t)] 



(13) 



This introduces extra integrations over the virtual momentum q increasing considerably the numerical effort [7]. 
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1.2 




FIG. 1: (Color online) Relative probability PEx{h,t) / PE\{h,t — oo), for production of a muonic atom in pp collisions at the 
LHC with an impact parameter b — amuon ~ 255 fm, where flmuon is the Bohr radius for a muonic atom. 



Recently, it was shown that one can treat the effects of close collisions {R' < R) by using the potentials of eqs. 
(3][5| for R' > R and the potentials for close collisions given by [TOI (we will only treat the El case for reasons to be 
shown later) 



(r, b, t) = (r) (1 - l,l.)i,Y*, (f ) (r) Ze 



1 



{<?o(r) + c«52 (r)} 



y/2vt 



if K = 

if K = ±l ' 



where cq — 2/5, c±i ~ —1/5, and 

50 (r)- 4 



In [r(i + v)] , 32 (r) 



4ti V v'^ 



1 In [r(i + w)] 



15 
2^' 



(14) 



(15) 



The integrals over r now have to be carried out by using eq. ([3| for R! > R and eq. (14) for R! < R. In particular, 
for 5 = one gets 



dt'e'-*Vj\r (r, b, t') = <5,,o*« (r) (1 - lol.)^Y*^ (?) *p (r) W ^^e^Fz; {50 (r) + C052 (r)} 



sm 



Rlu\ Rio 
Tv J Tv 

Ruj 



cos 



Ruj\ 
I\7 j 



if i? < Tvt 



(16) 



Tv 



+ ie'*'^ {itLU ~ 1) if i?>rwi 



For F ^ 1, (7o(F) + cog2(r) = 21n(2F) — 15/2. When t — 00 only the upper term in the last part of this equation 
contributes {R < 00) to the production probability. 

We use the above separation of close and distant collisions to make an estimate of their relative contributions to 
the total cross section. To obtain the total cross sections we use the formalism developed in ref. [7] which allows a 
better reduction of the integrations in momentum space. Adapting their results for ion-ion collisions, we get 



cr = 



I de \M{e)\ - dq , , . , . , , , ^ • (17) 



Notice that we keep the relative velocity v between the colliding nuclei in all formulas, although v ^ c. This 
is necessary because sometimes important combinations of 1 and v conspire and combine into Lorentz factors F = 
(1 — ti^)~^/^ in subsequent steps of the calculations. 



III. RESULTS 



At the Large Hadron Collider (LHC) at CERN the Lorentz gamma factor in the laboratory frame, Tiat, is 7000 for 
p-p, 3000 for Pb-Pb collisions. The relationship between the Lorentz contraction factor associated with the relative 
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along beam axis 



FIG. 2: (Color online) Coordinate distribution in the plane along the beam axis of the probability for production of muonic 
atoms with capture in the K-shell with pp collisions at the LHC. The impact parameter is chosen as 6 = amuon = 255 fm and 
is along the vertical axis. 



velocity between the colliding nuclei, and the collider energy per nucleon, E/A, in GeV, is given by F = 2(F^^(j — 1) ~ 
2{1.0735E/A)^ . This means that for the production of exotic atoms we have effectively F ~ 10^ for p-p and F ~ 10^ 
for Pb-Pb collisions. 

In figure [l] we plot the relative probability Psxib, t)/PEi{b, t = oo), for production of a muonic atom in pp collisions 
at the LHC with an impact parameter h = amuon = 255 fm, where amuon is the Bohr radius for a muonic atom. The 
time scale is in units of 10~^ fm/c. The probability is normalized so that it is equal to the asymptotic value of the 
El multipolarity. In absolute values the probabilities are very small, justifying the use of perturbation theory. The 
E2 probability tends to its asymptotic value faster than El. This can be understood from Eqs. ^ and Q as the 
time-dependence of the E2-field is determined by a higher inverse power of time. The asymptotic value of the E2 
probability is about a factor 4 smaller than the El case. As a function of the impact parameter, the E2 probability 
decreases with an additional 1/b^ dependence, as compared to El. This yields cross sections for pair-production with 
capture due to the E2 field being much smaller than that with El. The probabilities and cross sections are also much 
smaller in the case of the Ml multipolarity, due to the factor {pz/m) in Eq. [5j This is also substantiated by the 



approximation in Eq. ( 11 ), which has a factor \/ 1 — m instead of the -^/e + m which appears in the electric multipole 
cases. For low positive particle momentum this leads to a further suppression of the Ml multipolarity. 

We thus conclude that the El multipolarity alone is responsible for most part of the bound-free pair-production 
probability. This comes as no surprise because of the very large F factor in the frame of reference of the exotic 
atom. The spatial distribution of the time-dependent field is compressed as a pancake-like object with a spatial 
width Az — cAt — b/T. For b — 200 fm this is equal to Az ~ 10~^ fm for the LHC. Even for very large impact 
parameters, e.g., lA = 10^ fm, Az is small compared to the nuclear sizes. The E2 field is a measure of the "tidal" 
force, proportional to the spatial spreading of the electric field [3]. This "tidal" effect becomes larger at smaller 
velocities, when the field lines are not as compressed. The large value of F also leads to a complete dominance of the 
K = ±1 component of the El field in eq. 

The above discussion implies that in order to calculate pair-production with capture in ultra-peripheral collisions 
at the LHC once only needs to consider the El field, with k = ±1, in eq. The asymptotic production probability 
amplitude, Eq. ([2]), for a given impact parameter 6, is then given by 

a.=±i(p,^) = 2«z-^yyCXi(0 j dV¥„(r)(l-7o7.)rri„(f)vI'p(r), (18) 

where ^ = lob/Tv and Ki is the modified Bessel function of first order. Note that £,Ki{S) ~ 1 for ^ < 1. For 
^ > 1, £,Ki(S^) drops to zero exponentially. This means that the production probability drops as 1/6^ up to bmax ~ 
Thc/m K, 10^°(MeV.fm)/mc^ for pp collisions at the LHC. For production of muons and pions with capture this 
means bmax ~ 10^ fm, whereas for proton-antiproton with capture, this means bmax ~ 10^ fm. 

In Fig. [2] we show the coordinate space distribution for production of muons with capture in the K-shell in pp 
collisions at the LHC. The distribution is shown in a plane containing the beam direction and the impact parameter 
vector. The calculation is done for the El multipolarity and for an impact parameter b = amuon- The positive muon 
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FIG. 3: (Color online) Same as the previous figure, but in the plane perpendicular to the beam axis. The impact parameter 
vector lies along the x direction. 
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FIG. 4: (Color online) Probability of muonic atom production in pp collisions at the LHC as a function of the impact parameter. 



energy is taken as e = l.lm^ and its direction of emission is chosen as 100° when measured along the direction of 
motion of the muonic atom. One notices that the production mechanism is more efficient in regions perpendicular to 
the beam axis. The darker areas are representative of larger production rates in coordinate space. Figure |3] shows 
the same result, but as seen in the plane perpendicular to the beam axis. In this case, the impact parameter vector 
lies along the x-direction. We observe that the probability density is largest within a torus-like region with a radius 
r « amuon from the origin of the muonic atom, with the torus axis along the beam direction. 

We now use Eq. ( 16 1 to obtain the production probability of muonic atoms at the LHC at zero impact parameter. 
We find that the probability is 5 orders of magnitude smaller than that with b = amuon- We easily understand this 
result by inspection of the equation (16 1. When t = oo only the upper term in the last part contributes (i? < oo) 
to the production probability. This term oscillates harmonically as a function of loR/Tv ~ 2mR/Tv. For R along 
the z-direction this variable varies as 2mc^z/hc. As the largest contribution to the integral arises for z k, amuom this 
variable is much larger than the unity, causing the harmonic functions to oscillate wildly, and leading to a small value 
of the integral over coordinates (matrix element). 



We now perform an estimate of the contribution of collisions with small impact parameters by using Eq. ( 16 ) and 



integrating it numerically over time, space, and impact parameters, for several energies and angles of the emitted 
positive muon. To obtain these estimates, we have used the functions ^ - and ^+ without the derivative corrections. 
The integrals are performed for b < amuon- We compare with the results obtained using Eq. (18) for b > amuon- We 
find that the impact parameter region b > amuon contributes at least a factor 100 more to the cross section than the 
region with b < amuon- Such findings are confirmed by an integration over the free positive muon energy. This is 
also confirmed using Eqs. ( 18 ) and ( 16 ) and and ^4. without derivative corrections. Our numerical results are 
shown in figure [4] We clearly see that the small impact parameters, b < amuon yield a much smaller probability than 
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FIG. 5: (Color online) Contour plot with the angular distribution of the positive muon when the negative muon is captured by 
a proton at the LHC, as a function of the angle that the free muon has with the direction of motion of the muonic atom and 
of the energy of the free positive muon. 



exotic atom 


PP 




Pb-Pb 


anti-hydrogen 


63.4 pb 




132 b 


muonic 


44.8 X 10 


pb 


0.16 mb 


pionic 


21.3 X 10 


pb 


0.09 mb 


kaonic 


1.3 X 10" 


4 pb 


4.3 ^lh 


p-atom 


0.51 X 10 


pb 


1.3 ^lh 


protonium 


0.09 X 10 


pb 


0.3 iih 



TABLE I: Cross sections for production of exotic atoms in pp and Pb-Pb collisions at the CERN Large Hadron Collider (LHC). 



In Fig. [5] we plot the angular distribution of positive muon when the negative muon is captured by a proton at the 
LHC as a function of the angle that the free muon has with the direction of motion of the muonic atom and of the 
energy of the free positive muon. The impact parameter is chosen as 6 = CLmuon- The units in the contour plot are 
arbitrary, with the darker areas being the region of highest probability. The plot shows that the higher the electron 
energy is, the more backward peaked the distribution becomes. In the frame of reference of the atom, the angular 
distribution of the positive particles is backward peaked, along the beam axis, opposite in the direction of motion of 
the nucleus capturing the muon. The higher the positive muon energy is, the more backward peaked the distribution 
becomes. The most probable energy of the free muon is non-relativistic, i.e., e ~ m. For such low energy muons the 
opening angle for emission of the free muon is larger. In a collider, a Lorentz transformation of these results to the 
laboratory frame implies that all particles are seen along the beam direction, within an opening angle of order of l/F. 
That is, all positive particles are seen along the same direction as the beam. 



Using Eq. (17 1 we calculate the cross sections for production of exotic atoms in pp and Pb-Pb collisions at the 



LHC. Our results are shown in Table I. Whereas the cross sections are small for pp collisions, they are by no means 



negligible for Pb-Pb collisions. This is due to the factor in Eq. (17), although a reduction of this Z dependence 
arises from the distortion factor A/". 
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IV. CONCLUSIONS 



Because of the very strong electromagnetic fields of short duration, new and interesting physics arise at the LHC. 
We have studied the space-time dependence of the cross sections for production of exotic atoms in ultraperipheral 
collisions at the LHC. We have considered the case of muonic, pionic and anti-protonic (or protonium) atoms. A 
very transparent and simple formulation is obtained using the lowest-order corrections for the positron and electron 
wavefunctions and a multipole expansion of the electromagnetic field separated in the regions of small (i.e., b smaller 
than the Bohr radius) and of large impact parameters. Whereas most of our discussion and conclusions have been 
based on calculations of muonic atoms, the qualitative aspects will not change for pionic, protonic, or other exotic 
atoms, except for the magnitude of the cross sections. The space-time and impact parameter dependences will also 
be similar for Pb-Pb collisions. 

We considered the capture of negative particles in the K-shell of either the proton of the Pb nucleus. This is the 
largest contribution to the capture cross section. Inclusion of capture to all other shells will increase by about 20% of 
the our calculated cross sections, according to early estimates for electron-pair production with capture [3j. Depending 
on the relevance of this process to peripheral collisions with relativistic heavy ions and on the accuracy attained by 
the experiments, further theoretical studies on the capture to higher orbits could be necessary. The distortion of the 
wavefunctions due to the nuclear size should also be considered. Due to their small charges, pp collisions will not yield 
a measurable amount of exotic atoms, but in Pb-Pb collisions we expect an abundant number of events of production 
of exotic atoms. This will open opportunities to study the properties of exotic atoms and their decay widths. 
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